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Abstract 



The unintegrated quark and gluon distributions in the pion and nucleon are analysed 
using the CCFM equations in the single-loop approximation. We utilise the transverse- 
coordinate (or impact-parameter) representation which diagonalises the equations and 
study in detail the dependence on the transverse coordinate induced by the CCFM evo- 
lution. We find considerable broadening of the transverse-momentum distributions with 
an increasing magnitude of the hard scale, Q. For instance, at the Bjorken x = 0.1 the 
root mean squared transverse momentum of the gluons is enhanced by about 1 GeV when 
evolved from the the initial low scale to Q 2 = 10 GeV 2 , and by about 2 GeV when evolved 
up to to Q 2 = 100 GeV 2 . The broadening effect is enhanced with decreasing Bjorken x, 
and is stronger for the gluons than for the quarks. Analytic solution for the average 
transverse momentum corresponding to the x-moments of the distributions is obtained. 
The parton luminosities are also discussed. 



1 Introduction 



Parton distributions are fundamental and universal quantities of the QCD-improved parton 
model. They depend on the longitudinal momentum fraction of the parent hadron, x, as well 
as on the hard scale of the process, Q. The inclusive cross-sections describing hard processes are 
governed by the (collinear) factorisation theorem expressing those cross-sections as convolutions 
of partonic distributions and hard (partonic) cross-sections pQ . On the other hand, less inclusive 



measurements are often sensitive to the transverse momenta of the partons and require an 
introduction of the unintegrated parton distributions, i.e. distributions which are not integrated 
over the transverse momentum of the parton El HI El El El El El CHI HI]- These unintegrated 
distributions depend then upon the transverse momentum of the parton besides the dependence 
on x and Q. Typical cases where the information on the transverse momentum of partons is 
needed are the prompt photon production, or the transverse momentum distributions of the 
Drell-Yan lepton pairs, and Z° bosons, etc. The unintegrated parton distributions are also 
needed in the calculations made within the kx- factorisation framework [T2*|ITBtn"4U15| ll6j . 

The unintegrated parton distributions are described by the Catani-Ciafaloni-Fiorani- Mar- 
chesini (CCFM) equations [TT|ll8|ll9[l20p2ip221l23p24p25|l26|l27| which are based on the quantum 
coherence implying angular ordering along the partonic cascade [2E]- in the so-called single- loop 
approximation [T$ll2Ti] . adequate for large and moderately small values of x (say, x > 0.01), the 
CCFM equations reproduce the conventional leading-order DGLAP equations for the integrated 
distributions. Thus, they may be viewed as an extension of the DGLAP evolution which allows 
to investigate more general quantities. 

The purpose of this paper is to analyse in detail the unintegrated parton distributions of 
pions and nucleons. To this aim we explore the system of the CCFM equations in the single 
loop approximation utilising the transverse-coordinate representation of the parton distribu- 
tions. The transverse-coordinate, denoted in this paper by b, is a variable conjugate to the 
transverse momentum of the parton. It has been widely used in the study of the soft gluon 
resummation effects in the e + e~ collisions [2lHl5D| . in the transverse-momentum distributions of 
the Drell-Yan pairs or gauge bosons [31] , etc. The formalism adopted in our analysis is similar 
to methods used in those studies. An important merit of our use of the transverse-coordinate 
representation is the fact that it diagonalises the system of the CCFM equations in the single- 
loop approximation [321EIII- In practice, this means that the equations are solved independently 
for each value of h. Also, in this approximation the non-perturbative effects encoded in the 
choice of the initial dependence on b at the reference scale Qq (the profile function) factorise 
in the solution at a scale Q. This input profile introduces the damping factor for large values 
of b that makes it possible to extend smoothly the partonic /^-distributions down to the point 
fiT = 0. The CCFM evolution changes the dependence on b and leads to significant broadening 
of the transverse momentum distributions. The average transverse momentum squared grows 
approximately linearly with Q 2 (modulo logarithmic corrections). The effect is strongest for 
gluons, when e.g. for Q 2 = 100 GeV 2 the CCFM evolution is found to generate (k^Uli about 
2.3 GeV for x — 0.1 and 3.2 GeV for x = 0.01. These numbers should be compared to the initial 
spreading, (k^)p, which is typically assumed to be in the range of (0.5 — 1 GeV). Broadening 
of the k,T distribution corresponds to the development of the long-range tail in the distribution 
at large k T . 

The average transverse momenta squared can be expressed in terms of the logarithmic 
derivative of the distributions in the b representation at b = 0. Factorisation of the input 
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profile implies that the average transverse momentum squared of the parton is the sum of the 
'primordial' (non-perturbative) component at Qo and the term generated by the evolution. 
The unintegrated parton distributions are constrained by the integrated distributions which 
correspond to the integrals of the unintegrated distributions over the transverse momentum 
squared of the partons. In order to make our analysis realistic we use constraints implied by 
the results of the global LO QCD analysis jSUESl for the integrated distributions of pions 
and nucleons. To be precise, we take the same parametrisation of the starting integrated 
distributions as those used in p£fll33] and assume, as an educated guess, a Gaussian transverse 
momentum distribution, uniform for gluons and quarks 1 . Combining this assumption with 
the evolution from the CCFM equations in the single-loop approximation one should get a set 
of realistic unintegrated parton distributions hadrons. These can be used in phenomenological 
applications in the mentioned region of large and moderately small values of x, where the 
single-loop approximation is expected to be adequate. 

The content of our paper is as follows: In the next section we recall the system of the CCFM 
equations in the single-loop approximation and discuss its transverse coordinate representation. 
In Sec. 3 we discuss the unintegrated parton distributions of pions and nucleons which follow 
from the numerical solution of the CCFM equations in the single-loop approximation. We 
present our results for the b profiles and for the transverse-momentum distributions, as well 
as discuss the increase of the average transverse momentum with the increasing magnitude of 
the hard scale, Q, at different values of x. In Sec. 4 we discuss the moments of the b profiles 
and present results for the average transverse momenta squared, < fcf,(n, Q 2 ) >j, for the nth 
moment for partons of species i. We also derive (semi) analytic expressions for < k^(n,Q 2 ) >j. 
In Sec. 5 we discuss partonic luminosities and, finally, in Sec. 6 we give a summary of our 
results. 



2 CCFM equations in the single loop approximation 



The original Catani-Ciafaloni-Fiorani-Marchesini (CCFM) equation ^7J^] for the uninte- 
grated, scale-dependent gluon distribution fc(x, k T , Q), which is generated by the sum of ladder 
diagrams with angular ordering along the chain, has the following form: 



f G (x,k T ,Q) = fy(x,k T ,Q) 



f ^6(Q - qz)Q(q- q )^-A s (Q,q,z) 



ixq« J x z 

2N z fx 

2N c A NS (k T , q, z) + n \J {-, \kt + (1 - z)q\,q 



(1-z) 



(1) 



x The assumption of uniformity may be lifted with no difficulty, as well as other initial parameterisations 
than |.'llU.'i5| may be tested. This has little influence on our general conclusion on the large dynamical spreading 
of the distributions. 
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where A S (Q, q, z) and A NS (kT, q, z) are the Sudakov and non-Sudakov form factors, 



A s (Q,q,z) 
A NS (k T ,q, z) 



cxp 



cxp 



Q2 dp^Os_ /•' ' 

{qz)2 p 2 2ix J 
1 dz' f k r dp 2 2N c a s 



dzzP gg {z) 



27T 



(2) 
(3) 



The variables x, kr, and Q denote the longitudinal momentum fraction, the transverse mo- 
mentum of the gluon, and the hard scale, respectively. The latter is defined in terms of the 
maximum emission angle [TT)|I17| . The constraint 0(Q — qz) in Eq. (0) reflects the angular order- 
ing, and the inhomogeneous term, f°(x,k T ,Q), is related to the input non-perturbative gluon 
distribution. It also contains effects of both the Sudakov and non-Sudakov form-factors 



Equation Q in a sense interpolates between a fragment of the DGLAP evolution at large x 
and the BFKL dynamics at small x. To be precise, it contains only the g —>■ gg splittings and 
only those parts of the splitting function P gg (z) which are singular at either z —>■ 1 or z — > 0. 
In order to make the CCFM framework more realistic in the region of large and moderately 
small values of x we make use of the following extension [SS] : 



1. We introduce, besides the unintegrated gluon distributions fa{x, kr, Q), also the uninte- 
grated quark and antiquark gluon distributions, f qi (x, k T , Q) and fq^x, kx, Q)- 

2. We include, in addition to the g — ► gg splittings, also the the q — > gq, q — ► gq, and g — > qq 
transitions along the chain. 

3. We take into account the complete splitting functions P a j,(z), and not only their singular 
parts. 

In the region of large and moderately small values of the parameter x (x > 0.01 or so) one can 
introduce the single-loop approximation that corresponds to the following replacements: 

Q(Q-qz) - Q(Q-q), 
A NS {k T ,q,z) -> 1. (4) 

It is also useful to 'unfold' the Sudakov form factor such that the real emission and virtual terms 
appear on an equal footing in the corresponding evolution equations. Finally, the unfolded 
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system of the CCFM equations in the single-loop approximation has the following form: 



f l NS (x,k T ,Q) = f™ s (x,k T ) + / dz 



d 2 q a s (q 2 
7rq 2 2tt 



(5) 



x 



®( z - x )Pns \ -i k T,q - Pns( x , k T, q) 



fs(x,k T ,Q) = f° s (x,k T )+ Cdz f ^l^le(q 2 - q 2 )Q(Q - q) 

Jo J nq z 27T 



(6) 



x < Q(z — x) 



X 



P qq {z)fs \ -,k' T ,q) +P qg (z)f G (-,k' T ,q) - P qq (z)f s (x,hr,q) 



X 



f G (x,k T ,Q) = f° G ^k T ) + j\zj^^-e(q 2 -q 2 )e(Q-q) 



(7) 



x \Q{z-x) P gq (z)f s (-,Q' t ,q)+P gg (z)f G l-,Q' t ,q 



x 



(z) + zP qg (z)]f G {x, kr,q) \, 



where 



, T _ ^ T -r (1 - z)q. (8) 

The functions f* NS (x,kT,Q) are the unintegrated non-singlet quark distributions, while the 
unintegrated singlet distribution, fs(x, kx,Q), is defined as 

f s (x, k T , Q) = Z{ =1 [f qi (x, k T , Q) + f qi (x, k T , Q)] (9) 

The functions P a b(z) are the LO splitting functions corresponding to real emissions, i.e.: 

Al + z 2 



3 1-z 

2 



P qg (z) = N f [z 2 + (1 - z)\ 

4 1 + (1-z) 2 

3 ~z ' 

z 1 — z 



P 9q (z 



2N r 



1-z 



+ 



+ z(l-z) 



(10) 

(11) 
(12) 

(13) 



where Nf and N c denote the number of flavours and colours, respectively. After integrating 
over d 2 k T on both sides of Eqs. (0 - EI) we get the usual DGLAP equations for the integrated 
parton distributions pk(x,Q 2 ), defined by 

poo 

(14) 



(15) 
(16) 



xp k (x,Q 2 )= dk^f k (x,k T ,Q). 
Jo 

Equations ©-© can be diagonalised by the Fourier-Bessel transform, 

POO 

f k {x,k T ,Q) = / dbbJ (k T b)f k (x,b,Q), 
Jo 

POO 

f k (x,b,Q) = / dk T k T J (k T b)f k (x } k T ,Q), 



where k = NS, S, or g, and Jq(u) is the Bessel function. At b = the functions f k (x, b, Q) are 
related to the integrated distributions Pi(x, Q 2 ), 



f k (x, b = 0, Q) = -xp k (x, Q 2 ). 



(17) 



The corresponding evolution equations for f^s(x,b,Q), fs(x,b,Q), and f g (x,b,Q), equivalent 
to Eqs. © - © read 



2 dfNs(x,b,Q) 



Q 



2tt 



3Q 2 



dzP qq (z) 



x 



e(z-x)J [(l-z)Qb]f NS (-,b,Q)-zP qq (z)f NS (x,b,Q) , (18) 



2 df s (x,b,Q) 



Q 



2tt 



<9Q 2 



dz{ Q(z - x) J [(l - z)Qb\ P qq (z)f s i-,b,Q)+ P qg (z)f G l-,b,Q 



-[zP qq (z) + zP gq (z)}f s (x,b,Q) 
Q 



(19) 



2 df G (x,b,Q) 



a s {Q* 



2tt 



dzle{z - x) J [(l - z)Q6] [> OT (z)/s (|, 6, Q) + P ot (z)/g , 6, Q 

-[^P gg (2;)+zP gg (z)]/ G (x,fo,g) 
with the initial conditions 

f k (x,b,Q ) = f° k (x,b). 



(20) 



(21) 



In our analysis of the CCFM equations we assume, for simplicity and from the lack of 
detailed experimental knowledge, a factorisable form of the initial conditions ()21j) 



f° k (x,b) = -F(b)x Pk (x,Ql). 



(22) 



Other, more complicated and non-factorisable initial conditions may be explored with no diffi- 
culty as well. The input profile function, F(b), is linked through the Fourier-Bessel transform 
(llfjjl to the non-perturbative kx distribution at the scale Qo- At b = we have the normalisation 
condition for the profile function, F(0) = 1. 



3 Unintegrated parton distributions in the pion and nu- 
cleoli 

Since the initial condition for the QCD evolution, Eq. (|22p . assumes a uniform dependence 
on the profile function F(b) for all distribution functions fc, it is convenient to factorise this 
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function by introducing the 'tilde' distributions, 

fi(x, b, Q) = fi(x, b, Q)/F(b), z = G, S, NS. (23) 

The evolution equations (|18H2(J|) can then be written in terms of the / functions only. In other 
words, all dynamical information on the 6-dependence, generated by the evolution and linked 
to the Q and x variables, is contained in the / functions, while the uniform initial profile F(b) is 
carried over as a multiplicative factor and is decoupled from the dynamics. Certainly, we may 
choose different profiles F(b), which affects directly the shape of the unintegrated distribution 
functions, nevertheless the evolution is independent of this choice, as long as the factorisation 
is assumed. For this reason we first present the numerical results for the scaled functions 
fi(x,b,Q). 

The evolution equations for the / functions are solved numerically. The method used is 
based on the discretisation made with the help of the Chebyshev polynomials and is an exten- 
sion of the method developed in (HE] for the solution of the DGLAP equations. 

Figure ^ shows the solutions of Eqs. (jl8H20J) for the case of the pion with the LO Gliick- 
Reya-Schienbein (GRS) distributions jMj taken at the initial scale Ql = 0.26 GeV 2 . The left 
and right sides of the figure display the results for x = 0.1 and x = 0.01, i.e. moderately 
large and moderately small values of the Bjorken variable. We show, from top to bottom, 
the distributions fa, fs, and f^s, plotted as functions of the transverse coordinate, b. The 
various types of lines indicate consecutive values of the momentum Q, with the solid line 
corresponding to the initial scale of the GRS parameterisation, Q 2 = Ql = 0.26 GeV 2 |34j, 
and the dashed, dash-dotted and dotted lines denoting results for Q 2 = 1 GeV 2 , 10 GeV 2 , and 
100 GeV 2 , respectively. The original distributions at Q = Qo are flat in b, since we use the 
rescaled functions ([23)1 . We note immediately from the figure that the increase of Q results 
in sharp peaking of the distributions at b = 0, with the effect stronger at the smaller value of 
x. The shrinkage of the distributions in b is most prominent for the gluons. For instance, for 
the gluons at x = 0.01 the curves drop to half of the values at b = around b = 2, 1.2, and 
0.6 GeV" 1 for Q 2 = 1, 10, and 100 GeV 2 , respectively. To the extent these distributions can be 
approximated with a Gaussian, that would correspond to widths in the transverse momentum 
of about 0.8 GeV, 1.4 GeV, and 3 GeV, respectively. The transverse-momentum distributions 
will be discussed in detail in the following parts of the paper. At the values of x displayed in 
Fig. Ethe shrinking effect in b is similar for the singlet and non-singlet quark distributions. 

We note that for the case of gluons at x — 0.01 the function fa becomes negative above 
b ~ 4 GeV" 1 . To our knowledge, this poses no physical problems. Positivity is required for 
the unintegrated distributions in the transverse-momentum space, since then the positivity of 
cross-sections, which are physical quantities, is guaranteed. The Fourier-Bessel transform of a 
positive function, however, need not be positive itself. We will see shortly that in the transverse- 
momentum space the distributions are positive. The effect of negative fo at larger b and small 
x is an artifact of the valence-like input for the gluons in the parameterisations of Refs. [HUES], 
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and can be understood as follows: At small x the evolution equation (J2UJ) is dominated on the 
right-hand side by the gluon contribution, hence can be approximated with the following form: 



2 df G (x,b,Q) 



Q 



a s (Q 2 
2tt 



dQ 2 



dzJ [(l - z)Qb\ {p gg (z)f G (-, b, Q) - [zP 99 (z) + zP qg (z)]f G (x, b, Q) 



(24) 



-{zP qg (z) + zP gg {z) [1 - J ((l - z)Qb)]f G (x, b, Q)} 

At small x the dominant region of integration in the first term on the right-hand side of Eq. 
is x < z « 1, hence we may approximate 1 — z by 1 in the argument of the Bessel function in 
this term. The solution of Eq. (J24|) with the initial condition 

1 



then reads 



f G (x,b,Q ) = -xg(x,Q ) 



f G (x,b,Q Q ) = ^T g (b,Q)xg{x,Q 2 cS {b,Q)), 



where 



T g (b,Q) = exp 




q2 dq 2 a s (q 2 



[ dz[zP qg (z) + zP gg (z)[l-J ((l 
Jo 



q 2 2tt 

and Ql s (b, Q) is obtained from the following implicit equation: 

< B (b,Q) dq 2 a s (q 2 ) f Q " dq 2 a s (q 2 ) 



(25) 
(26) 

z)qb)}} \ (27) 



q 2 2tt 



q 2 2tt 



Mqb). 



(28) 



The function xg(x, Q 2 ) is the integrated gluon distribution corresponding to the solution of 
equation (|2^j) at b = with the initial condition xg(x, Qq). At sufficiently large b the integral 
on the right-hand side of Eq. (|2~%j) becomes negative, which gives Ql s {b, Q) < Ql, and the 
effective 'backward' evolution of the valence gluons gives a negative xg(x, Q 2 S ) for sufficiently 
small x and for Ql s (b,Q) sufficiently smaller than Q^. This is the technical reason for the 
behaviour seen in Fig. ^ 

The mean squared transverse momentum for a given distribution is equal to 



(krp^X, Q ))i 



/ d 2 k T k 2 T fi{x, k T ,Q) 



-4 



d/db 2 f t (b,x,Q) 



f d 2 k T fi(x, k T ,Q) 

where the last equality follows from the expansion J (z) 
the decomposition ()23|) we may immediately write 



fi(b,x,Q) 



(kj>(x, Q ))i 
(k 2 T ) F 



(k^(x,Q 2 )) 

i,evol 



(k 2 T ) F + {k^(x,Q 2 ))^ emX) 
. d/db 2 F{b) 



F(b) 



b=0 



-4 



d/db 2 Ub,x,Q) 



f(b,x,Q) 



(29) 

6=0 

. Through the use of 
(30) 



(31) 



b=0 



which means that the total mean squared transverse momentum is the sum of the contribution 
from the profile F, which is constant and corresponds to the width at the initial scale Q = Qo, 
and the piece (kj)i >e voh which is due entirely to the evolution and is independent of the profile 
F . Because of these features we concentrate on {kj) ewo \ in the discussion below. 

In Fig. El we show (fcf.)^, plotted as a function of Q 2 , for x=0.1 (top) and x=0.01 (bot- 
tom). The solid, dashed, and dot-dashed lines denote gluons, singlet, and non-singlet quarks, 
respectively. Again, we clearly note the spreading effect with increasing Q 2 . The effect is 
strongest for gluons, which at Q 2 = 100 GeV 2 achieve {k^) l J^ ol about 2.3 GeV for x = 0.1 and 
3.2 GeV for x = 0.01. However, even at moderate Q 2 the spreading effect is sizeable, e.g. at 
Q 2 = 10 GeV 2 and x = 0.01 we find (A£)^ = 1.5 GeV, 1.0 GeV, and 0.8 GeV for gluons, 
singlet, and non-singlet quarks, respectively. The quoted numbers should be compared to the 
initial spreading, (k^)p, which is typically assumed to be in the range of (0.5 — 1 GeV) 2 . We 
also note from Fig. (J2J) that at larger values of x the quark singlet and non-singlet widths are 
practically equal, and at lower x the singlet distribution is somewhat wider. 

We also note a very sharp growth of (fex)g^ ol at low Q and low values of x. This feature 
is an artifact of the valence-like gluon distributions at the reference scale Qq. 

In Fig. El we show the effect of the CCFM evolution on the kx distributions. We plot the 
unintegrated distributions fa,s,Ns( x , &t, Q) in the pion for x = 0.1 and x = 0.01, and for four 
values of the hard scale Q 2 varying from Qq to Q 2 = lOOGeV 2 . We have assumed here the 
Gaussian profile, 

F(b) = exp(-b 2 k 2 /A), (32) 

which leads to Gaussian kr distributions at the the initial scale, Qo, proportional to exp(— k^/k^). 
The width parameter is taken to be k^ = lGeV 2 . The broadening of the distributions with 
increasing magnitude of Q 2 and decreasing x is clearly visible. It can also be seen that the effect 
is strongest for the gluons. Important property of the broadening is significant modification of 
the Gaussian shape and development of the non-Gaussian long-range tail at large k T . 

Strong increase of the average transverse momentum squared with increasing Q 2 does also 
indirectly affect Q 2 dependence of unintegrated distributions at k T = 0. In order to understand 
this dependence one may use the following approximate relation: 



*(*.fr = o,q)~ < ggff ) ' >< (33) 



where Pi(x,Q 2 ) are the integrated distributions. (Equation (j^Hj) follows from approximating 
in (JT5j) the distributions fi(x,b,Q) by a Gaussian in b 2 .) One expects that fi(x,kr = 0,Q) 
should decrease at large Q 2 due to strong increase of < k^x^Q 2 ) >i with Q 2 that can be 
clearly seen from the from Figure particularly for the gluons. In fact we show in the next 
Section that < k^(x, Q 2 ) >« should increase as Q 2 (modulo logarithmic effects). This increase is 
stronger than the potential increase of xpi(x, Q 2 ) caused by scaling violations. Initial increase of 
fc(x, kx = 0, Q) with increasing Q 2 for low values of Q 2 , which can be clearly seen, particularly 
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for low value of x — 0.01 is caused by the fact that at low Q 2 < k^(x,Q 2 ) > eV oi,G« k 2 ,, i.e. 
< k^(x, Q 2 ) >g~ k\ and so the Q 2 dependence of fc(x, kr = 0, Q) is entirely driven by scaling 
violations of xp g (x, Q 2 ) = xg(x, Q 2 ). 

Figures EJIHl show the results analogous to Figs. for the case of the nucleon, where the 
Gliick-Reya-Vogt parameterisation of Ref. [US] is used at the initial scale of Ql = 0.26 GeV 2 . 
The results are very similar to the case of the pion presented above. For the case of the 
transverse-coordinate distributions we provide in Fig. 01 the results for both the u and d valence 
quarks. For the nonsinglet distributions the other figures show the results for the u valence 
quarks only, since the case of the d valence quarks is very similar. 



4 Evolution of the x-moments 

The evolution equation for the x-moments of the non-singlet distribution, 



fas(n,b,Q 2 )= [ dxx n - l f m {x,b,Q), (34) 
Jo 



is particularly simple. It immediately follows from Eq. (JTSjl that 

Q , dUnAQ 2 ) = as^l £ dzPqq{z) [z nj Q ((1 _ z)Qh)) _ !] /Ns(n , b} Q)} 



(35) 



which can be transformed into 



aiog(/ NS (n,6,<3 2 )) a s (Q 2 ) [' . 4z 2 + l 



V h*? = ^ I d ~~— ((1 " '' )Qb)) ~ 11 ' (36) 

At b = this equation reproduces the standard evolution of non-singlet moments for the 
DGLAP equations. Next, we differentiate with respect to b 2 at b = on both sides of Eq. (}3T)j) 
and evaluate the z integral on the right-hand side. The simple result is 

<9(/4(n,Q 2 )) NS ,evoi _ a s (Q 2 ) , , 



with the definition 



and 



dQ 2 2tt 



(kl(nm^= S ^tf^ ] (38) 

J d 2 k T f NS (n, k T ,Q 



4/1 1 1 1 , 

Cn ~3\n + l~n + 2 + n + 3~n + 4'° ()<)i 
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The integration of Eq. (J37)l with the condition (kj,(n, Qo))NS,evoi(<5o) = yields the final result 



f Q2 a (a 2 ) 

{k 2 T (n, Q 2 )>NS,evoi = c n / dq 2 ^l, (40) 



where 




with /3q = lliVc/9 — 2Nf/3 and li(^) = J* dt/logt. Obviously, the moments for all values of n 
are growing with Q 2 , which once again is a manifestation of the advocated spreading in k^. 

The analysis for the gluon and singlet quarks is more complicated due to mixing and is 
discussed in the Appendix. The evolution parts < k^(n, Q 2 ) >i iev oi, with i — G, S, are equal to 

„ l2f ,o2\ ^ A fi,cm\{ n iQ ) /.(-.n. 

< k T {n,Q ) >i,cvoi= -4- — (42) 

/i,evol(n, Q Z ) 

with the functions fl evo i(n, Q 2 ) and fi !e voi(n, Q 2 ) defined in Eqs. (f5%|) and (f72|) in the Appendix. 
From those formulas and from equation (|4T)|) it follows that < k^(n,Q 2 ) >; increases as Q 2 , 
modulo logarithmic effects. This increase is caused by the fact that —f' evo \(n,Q 2 ) is found to 
increase as Q 2 modulo logarithmic effects caused by QCD scaling violations. Similar strong in- 
crease should also hold for —d/db 2 fi(b, x, Q)\b=o controlling Q 2 dependence of < k\{x, Q 2 ) >i tevo i 
The results for the dynamically-generated root mean squared transverse momentum for the x 
moments of the distributions in the pion, < fc^(n, Q 2 ) >i jev oi, are shown in Fig.[7| For n = only 
the nonsinglet distribution is shown, since the gluon and singlet equations involve singularities 
and do not make sense in this case. The solid, dashed, and dot-dashed lines denote gluons, 
singlet quarks, and non-singlet quarks. Again we use the initial GRS condition 34 j at the scale 
Q 2 = Ql = 0.26 GeV 2 . The analogous plot for the nucleon with the GRV parameterisation [HH] 
gives identical results for the nonsinglet distributions, which follows immediately from Eq. (JHHj) 
not involving the initial values of the moments. For the gluons and singlet quarks the results 
for the pion and the nucleon, although not strictly equal (see the Appendix), are practically 
identical. The results displayed in Fig. [7| once again show the main observation of this paper, 
namely, a large spreading of the distributions in process of the single-loop CCFM evolution. 



5 Parton-parton luminosity function 

The knowledge of unintegrated distributions in the b representation is useful for calculations 
of the parton-parton luminosity function corresponding to the collision of parton a with the 
longitudinal momentum fraction x\ and the transverse momentum fraction k\T with parton b 
with the longitudinal momentum fraction X2 and the transverse momentum fraction /c2T- The 
luminosity function is defined as 0: 

L ab (x u x 2 ,q T iQ) = / — —fa(x 1 ,k 1T ,Q)f b (x 2 ,k 2 T,Q)5 i2) (kiT + k 2T - q-r)- (43) 
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Substituting in Eq. the following representation of the 5 function, 



IT + k2T — qTj 



(2vr) 2 



1 



/ 



<i 2 bexp[(ki T + k 2 T — qiOb] 



(44) 



we find that the integral defining the parton-parton luminosity in Eqn. (J43)) can be directly 
expressed in terms of the unintegrated distributions in the b representation, 



In Fig. |H1 we show the luminosity function corresponding to the valence u quark-gluon colli- 
sions. This quantity is relevant for the description of prompt photon production for which the 
dominant subprocess is u + g — » u + 7. We set Q = p 7 x and x\ = x 2 = Zp^x/ \fs where \fs is 
the CM energy of the pp collisions. The two curves correspond to y/s = 24.3GeV, p^x = 4GeV 
and \fs = 1800GeV, p 7 x = lOGeV. We can see that the luminosity function is increasing 
with the energy, which simply reflects the increase of the unintegrated gluon distributions with 
decreasing x = 2p lT / \fs. 

6 Summary and conclusions. 

In this paper we have performed an analysis of the unintegrated parton distributions fi(x, kx, Q) 
of the pion and nucleon using the CCFM equation in the single loop approximation. For the 
integrated distributions this approximation is equivalent to the LO DGLAP evolution, and thus 
should be adequate in the region of large and moderately small values of x {i.e. x > 0.01 or 
so). Important merit of the CCFM framework is the fact that once the input x and kx parton 
distributions are provided at the reference scale Ql, the CCFM evolution generates the x and 
hr distributions for arbitrary Q 2 > Qf, and arbitrary values of kx, including the region of low 
kx down to the point k^ = 0. This framework does therefore provide a useful extension of the 
approximate DGLAP - based treatment of the unintegrated distributions discussed in - 
We have extended the original CCFM equation by including the quarks and non-singular parts 
of the splitting functions and used the transverse-coordinate representation by introducing 
the parton distribution fi(x, b, Q), related to the unintegrated parton distributions fi(x, kx, Q) 
through the Fourier-Bessel transform. The usefulness of this representation is related to the 
fact that the transverse coordinate b is conserved through the CCFM evolution in the single 
loop approximation. The average transverse momenta squared of the partons are determined 
by dfi(x,b,Q)/db 2 \ b=Q . We have studied the impact of the QCD evolution on the shape of 
the b profiles and on the broadening of the kx distributions. We have quantified increase 
of the average transverse momentum with increasing magnitude of the hard scale Q and/or 
with decreasing x and we gave semianalytic insight into this increase. The average transverse 
momenta squared of the partons were found to increase as Q 2 with increasing Q 2 , modulo 
logarithmic effects related to the QCD evolution. We have also pointed out that the parton 




(45) 
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distributions in b representation determine the parton luminosity functions for fixed sum of the 
transverse momenta of the partons. The unintegrated distributions fi(x, kx, Q) (or fi(x, b, Q)) 
obtained in our paper can be used for the theoretical description of the processes which are 
sensitive to the transverse momenta of the partons. 
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Appendix 



The average value of the transverse momentum < k^(n) >i is given by 

J>,Q 2 ) 



< k l T {n,Q 2 ) >i= -4 



Z*(n,Q 2 ) 



(46) 



where 



Ji{n,Q ) = ^2 l b2=0 ' 

= f t (n,b 2 = 0,Q 2 ), 



fi(n,Q 2 ) 

and fi(n,b 2 ,Q 2 ) are the moment functions, 

f l (n,b 2 ,Q 2 )= f dxx n ~ l UxAQ)- 
Jo 



(47) 



(48) 



The index i corresponds to the gluon, quark singlet, or quark non-singlet distributions. In the 
case of the quark NS distributions we just have 



f 1 f Q2 

<k^(n,Q 2 )> NS = <k 2 T {n,Ql)> m + dzz n (l - z) 2 P qq (z) dq' 

Jo Jol 



2ir 



= < k 2 T {n, Ql) > NS + < k 2 T {n, Q 2 ) > N s,ovoi 



(49) 



In what follows we derive analytic expression for the dynamical component < fcf,(n, Q 2 ) >i iev oi 
of < k\{n, Q 2 ) >j corresponding to the gluon and quark singlet distributions. To this aim we 
start from the system of the evolution equations: 



Q 



,df'(n,Q 2 ) a s (Q 2 



dQ 2 



2tt 



P(n)f(n,Q 2 )-^P'(n)f(n,Q 2 ) 



(50) 



obtained from Eqs. (...) by differentiation with respect to b 2 . The components of the vector 
f'(n, Q 2 ) are f' s (n, Q 2 ) and f' G {n, Q 2 ), while the components of the matrices P(n) are equal to 
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the moments of the splitting functions defining DGLAP evolution of the moments of singlet 
and gluon distributions, namely, 



P qq (n) 



P 99 (n) 



dzP qq (z)[z n - 1], 



dz{zP gg {z)[z n - x -1]- zP qg }, 

dzPqg t gq(z)z . 



The explicit expressions for the functions Pij(n) are 



P qq (n) 
P q9 (n) 



P 9q (n) 



P 9q {n) 



4 
3 

N f 

4 
3 



n + 1 n + 2 

L 2 2 

+ 



2S(n) 



n + 1 n + 2 n + 3 



2 1 

+ 



n n + 1 n + 2 



11 
12 



S(n) 



n n + 1 n + 2 n + 3 



where 



S(n) 



f/2 



1 - Z n 
1 - Z 



- 1 

^ jfe' 
k=l 



Similarly, the components of the matrices P'(n) are defined as 



P' ab {n) ■ 

with the explicit expressions given by 



dz{l - z) 2 P ab {z)z n , 



PM 



P' m {n) 



4 
3 

N f 

4 " 
3 

6 



1 



1 



n + 1 n + 2 n + 3 n + 4 
14 7 6 



+ 



n + 1 n + 2 n + 3 n + 4 n + 5 



6 7 
+ 



4 1 
+ 



n n + 1 n + 2 n + 3 n + 4 



3 '' 



3 5 
+ 



5 3 
+ 



1 



n n + 1 n + 2 n + 3 n + 4 n + 5 

The dynamical components < k\{n, Q 2 ) >i jCVO \ of < k\{n, Q 2 ) are given by Eq. 
fi{n,Q 2 ) = f[ evol (n,Q 2 ) corresponding to the solution of Eq. (J5Uj) with the initial 
f.(n, Ql) = 0. The solution of Eq. flUD with the initial condition f'(n, Q 2 = Qg) = 
formally written as 



(51) 



(52) 
(53) 

(54) 



fevol( n i 



Ql 



dq 2 ^!+p- exp{P(n) [£(Q 2 ) - £,{<?)]} P'{n)f{n, q 2 ) 



(55) 

with 
condition 
can be 



(56) 
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where 

This representation gives, after some algebra, the following explicit expressions for /s,evoi( n 5 Q 2 ) 
and /g,cvoi(«>Q 2 ) : 

f's,G,evol{ n i Q 2 ) — 

-^{exp^We^l^W^O^^+^W^-v^R,^)] 

+ exp[ 7 -(n)e(Q 2 )] [^,g(^)^(0, Q 2 ) + ^(^^(a/AH, Q 2 )]}, (58) 

where 

^(r,Q 2 )= T rfg 2 ^expK(g 2 )] (59) 
The functions 7 mp (n) are the eigenvalues of of the matrix P(n) and are given by 



7*(n) = W + W^^ (60) 

with 

A(n) = [P w (n) - P gg {n)f + 4P g »P OT (n). (61) 
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The coefficients hg^,(n), hg G (n),h SG (n), and h s ^(n) are expressed in terms of the input 
moments fs,c(n, Ql) 

h + s + {n) = ^^{[7 + (n)-^>)][^(n)^(n,Qg) + ^(n)/J(n,Qg)] 

+ P q9 (n)[^ q (n)fi(n,Q 2 ) + P^(n)f^n,Ql)]} (62) 

h + G + {n) = -^^{P gq (n)[P^ q (n)fUn,Ql)+PMf^(n,Ql)} 

+ h + (n) ~ P m {n)] [P gq (n)f+(n, Ql) + P' gg (n)f+(n, Ql)]} (63) 

ht(n) = ^={[7 + (n)-P 99 (n)][P;(n)/ s (n,^) + P;(n)/j(n,^)] 

+ P qg (n)[P gq (n)fs(n,Ql) + P^ g (n)f G (n,Ql)]} (64) 

h G ~(n) = ^={P OT (n)[^(n)/5(n,^)+^(n)/o(n,C^)] 

+ [7 + (n)-P 99 (n)][P; 9 (n)/ 5 (n,g^) + P;(n)/ G (n,g^)]} (65) 

+ ^ w (n)[^(n)/5(n,C^) + ii(n)/o(n,Qg)]} (66) 

h G -{n) = -^^{P 9q (n)[P^ q (n)fs(n,Ql) + P^n)fa(n,Ql)\ 

+ [ 7 -(n)-P g9 (n)][P; 9 (n)/ 5 -(n,g2) + P;(n)/ G (n,g2 )]} (67) 

= -^^{[7-W-P 95 (n)][P>)/ 5 + (n,g^) + P>)/+(n,gg)] 

+ p^WtP^Wz+^g^ + p^^z+^g^)]} (68) 

^G + H = --7^{^(n)[^(n)/r(n,^) + i^(n)^(n,Qg)] 

+ [7"(n) - P„(n)] [P;,(n)/ 5 + (n, Qg) + P' gg {n)f^n, g 2 ,)]}, (69) 



where 



/s Qo) = T^{h T (n)-P gg ]f s (n,Ql) + P qg (n)f G (n,Ql)}, (70) 
/^(n,gg) = T ^={[^(n)-P qq ]f G (n,Ql) + P gq (n)f s (n,Ql)}. (71) 



The functions fs,G(n, Q 2 ) are given by: 

fs, G (n,Q 2 ) = /^,gg)exp[ 7 + (n)e(g 2 )] +/c G (n,gg)exp[ 7 -(n)e(g 2 )] (72) 
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Figure 1: Solutions of the CCFM equations in the single- loop approximation, Eqs. (jl8H20j) . for 
the unintegrated partonic distribution functions in the pion. The gluon, fc, quark singlet, f$, 
and quark non-singlet distributions, /ns, are plotted as functions of the transverse coordinate, 
b. The left and right parts are for x = 0.1 and x = 0.01, respectively. The solid lines correspond 
to the initial scale of the GRS parameterisation, Q 2 = Ql = 0.26 GeV 2 while the dashed, 
dash-dotted and dotted lines correspond to Q 2 = 1 GeV 2 , 10 GeV 2 , and 100 GeV 2 , respectively. 
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Figure 2: The dynamically-generated root mean-squared transverse momentum of unintegrated 
partonic distribution functions in the pion, plotted as a function of Q 2 . The top and bottom 
parts correspond to x=0.1 and x=0.01, respectively. The solid, dashed, and dot-dashed lines 
denote gluons, singlet quarks, and non-singlet quarks. The CCFM equations in the single-loop 
approximation with the initial GRS condition [33] at the scale Q 2 = Ql = 0.26 GeV 2 are used 
for the evolution. 
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x = 0.1 x = 0.01 
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Figure 3: Transverse- momentum parton distributions in the pion following from the solution 
of the CCFM equation in the single loop approximation. The gluon, fa, quark singlet, fs, and 
quark non-singlet distributions, Jns, are plotted as functions of the transverse momentum, k?- 
The left and right parts are for x = 0.1 and x = 0.01, respectively. The solid lines correspond 
to the initial scale of the GRS parameterisation, Q 2 = Ql = 0.26 GeV 2 while the dashed, 
dash-dotted and dotted lines correspond to Q 2 = 1 GeV 2 , 10 GeV 2 , and 100 GeV 2 , respectively. 
Parton distributions /j are in GeV -2 and transverse momentum kx is in GeV. 
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Figure 4: Same as Fig. ^for the case of the nucleon, with the initial GRV parameterisation 
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Q 2 [GeV 2 ] 

Figure 5: Same as Fig. H]for the case of the nucleoli, with the initial GRV parameterisation [35J. 
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Figure 7: The dynamically-generated root mean squared transverse momentum of unintegrated 
partonic distribution functions in the pion corresponding to the x moments of the distributions, 
plotted as a function of Q 2 for the moment index n = (top), n = 1 (middle), and n = 2 
(bottom). The solid, dashed, and dot-dashed lines denote gluons, singlet quarks, and non- 
singlet quarks. For n = only the non-singlet distribution makes sense. The CCFM equations 
in the single-loop approximation with the initial GRS condition [33] at the scale Q 2 = Q\ = 
0.26 GeV 2 are used for the evolution. The analogous plot for the nucleon with the GRV 
parameterisation gives identical results for the nonsinglet distributions, and practically identical 
for the gluons and singlet quarks. 
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Figure 8: The parton luminosity function corresponding to the valence u quark - gluon collisions 
plotted as the function of the transverse momentum equal to the vector sum of the transverse 
momenta of the colliding partons. The value of the hard scale Q is set to the transverse 
momentum p lT of the produced photon. The two curves correspond to two values of the total 
pp CM energy y/s and different values of p lT '- = 24.3GeV, p lT = 4GeV (dashed line) and 
Vs = 1800GeV, p lT = lOGeV (full line). Results for ^/l = 24.3GeV were multiplied by 10. 
The parton luminosity function L is in GeV -2 and transverse momentum is in GeV. 
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